Enhancing the photomixing efficiency of optoelectronic devices in the terahertz regime 
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A method to reduce the transit time of majority of carriers in photomixers and photo detectors to 
< 1 ps is proposed. Enhanced optical fields associated with surface plasmon polaritons, coupled with 
velocity overshoot phenomenon results in net decrease of transit time of carriers. As an example, 
model calculations demonstrating > 280x (or --^2800 and 31.8 fiW at 1 and 5 THz respectively) 
improvement in THz power generation efficiency of a photomixer based on Low Temperature grown 
GaAs are presented. Due to minimal dependence on the carrier recombination time, it is anticipated 
that the proposed method paves the way for enhancing the speed and efficiency of photomixers and 
detectors covering UV to far infrared communications wavelengths (300 to 1600 nm). 



Recent experimental observation of extraordinary op- 
tical transmission through a normal metal film having 
subwavelength holes [l| has generated intense interest in 
harnessing the underlying physics for photonic applica- 
tions. It is widely believed that quasi two-dimensional 
electromagnetic excitations tightly bound to the metal 
surface known as Surface Plasmon Polaritons (SPPs) are 
responsible for the observed near field enhancement of 
the optical radiation. Exponential decay of field com- 
ponents arising from SPPs (penetrating ~100, '^lO nm 
in dielectric and metal respectively), make them highly 
attractive for miniature photonic circuits and devices[2|. 
Currently, several potential applications of SPPs are be- 
ing explored including wave guiding 3] and near-field 
microscopy 4J. However, application of SPPs to enhance 
the speed and efficiency of photo detectors or photomix- 
ers has not been attempted due to the complexity of 
underlying physics and optimal device designs involving 
complicated dielectric-metal structures. 

In this letter we propose a simple interdigitated metal 
structure embedded in a suitable semiconductor ma- 
terial to harness SPPs for improving the speed and 
efficiency of photomixers and detectors. Such a de- 
vice can be realized with the existing material growth 

FT, 

techniques |3|. As an example, we demonstrate via model 
calculations, significant (> 280 x) improvement in the 
efficiency of a THz photomixer based on Low Temper- 
ature grown GaAs (LT-GaAs). Microwave to THz ra- 
diation sources in the 0.1 to 10 THz are being exten- 
sively studied for their application in communications, 
medical imaging, and spectroscopyja [3, IM y> UM- For 
many applications- compact, narrow-linewidth, widely 
tunable continuous wave sources are necessary. In par- 
ticular, to be used as a local oscillator in communica- 
tions systems the THz source should produce stable out- 
put power > 10 /J,WJM |ll|. Among the various tech- 
niques being pursued [fl, 1^ lul ll^ ll^ electrical down 
conversion of optical microwave sources in a suitable 
photomixerja, |lj |l3 are appealing due to the easy 
availability of tunable, narrow-linewidth, and stable solid 
state lasers. However, until now the output power from 
these photomixers is limited to < 10 /iW in the cru- 
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FIG. 1: Schematic of FP cavity (a) and a layer of interdig- 
itated metal lines (b) modelled using FDTD and FD sim- 
ulators. The FP cavity is formed by 3 or 4 dielectric lay- 
ers constituting DBRi, A*' pairs of low/high refractive index 
dielectric layers constituting DBR2, and absorbing layer of 
thickness D. The active volume is terminated by 15 cell 
UPML layers in FDTD implementation. The lines parallel 
to the X axis (width 200 nm) are connected to either sig- 
nal transmission lines or planar antenna fabricated on top of 
the absorbing layer (not shown). Through layer vias electri- 
cally connect identical layers of interdigitated lines (b) placed 
at various depths Zi within the absorbing layer. The over- 
all dimensions of the absorbing volume (excluding DBRs and 
via leads) are given by La;, Ly, D whereas number of in- 
terdigitated metal layers, their vertical positions, (measured 
from the Air-DBRi interface), the metal line (finger) thick- 
ness, pitch, width, and lengths of -l-ve and -ve electrodes are 
given by i, Zi, d, p, w, h, I2 as shown, a; = 0, y — 
is located at the center of top half of the interdigitated pat- 
tern in (b). When the top and bottom interdigitated struc- 
tures in (b) are connected in series as shown with the center 
electrode grounded, the resulting structure can be excited by 
lasers forming two quasi-elliptical spots. Such a radiation pat- 
tern is achievable by choosing proper laser propagation mode 
(such as TEMoi mode). The spacing between the two halves 
is adjusted to match the beam pattern of the specific laser 
system being used. 



cial THz range. In a conventional Metal-Semiconductor- 
Metal (MSM) photomixer ^14J, optically thick ( 100 nm) 
interdigitated metal lines are fabricated on a high resis- 
tivity semiconductor with subpicosecond recombination 
time of carriers {TeJTh for electrons/holes respectively) 
and high breakdown field limit. So far, only a few materi- 
als such as LT-GaAsjftll4J and Fe implanted InGaAsp^ 



have been shown to meet these requirements. From the 
photomixer theory, when two lasers (wavelengths Ai, A2, 
'^o = (-^i+A2)/2, and powers Pi = P2 = -Fb/2) with their 
difference frequency / = c(Aj~^ ~ ^2^^) in the THz range 
are incident on a dc biased photomixer [lj|, the THz wave 
output power Pf is given by 



P, 



/ 



■^Rl^p 



[1 + i2Tr fiRL + Rs)Cry 



(1) 



where ip is the photocurrent generated in the external 
circuit, c is velocity of light in free space, C is the ca- 
pacitance, Rs is small internal resistance of the metal 
structure, and Rl is the load resistance in the 72-200 
range ja Ha . Il4| . To accelerate carriers generated deep in- 
side the semiconductor, high dc voltage (Vb « 40 V) is ap- 
plied across the electrodes. This results in fields quickly 
exceeding the breakdown limit near the electrodes lead- 
ing to device failure. As the carrier transport is transit 
time (ttr) limited, i. e., ttr > f~^, recombination time 
T <C /~^ is required to recombine the carriers that do 
not contribute to ip before the beat cycle reverses |l2lll4J . 
Even if the photomixer is placed in a suitable optical cav- 
ity, due to strong reflection from thick metal electrodes, 
no carriers are generated directly below the electrodes 
where ttr will be small. In addition, subwavelength fea- 
tures of the metal lines produce strong near field diffrac- 
tion patterns that are not taken into account in conven- 
tional designs H E [li| . 

In the present technique, we exploit SPPs for generat- 
ing more carriers close to subwavelength normal metal 
electrodes embedded in a photoconductor layer sand- 
wiched between two Distributed Bragg Reflectors (DBR) 
that form a Fabry-Perot (FP) cavity. When coupled with 
velocity overshoot phenomenon reported in many photo- 
conductor materials m E El EE H 113, Ur of ma- 
jority of carriers is reduced to < 1 ps, i.e., Ur < /~^. 
For efficient use of SPPs, it is desirable to have thin 
(~10 nm) normal metal lines with subwavelength fea- 
tures distributed throughout the active volume (Fig la) 
in a manner that would collectively enhance the optical 
field intensity in the vicinity of the metal lines. This 
must however be accomplished without significantly in- 
creasing C and Rs ■ We show that for the proposed struc- 
ture (Fig 1) optimized for Aq = 850 nm, Pf is ^2800 
and 31.8 /xW at / = 1 and 5 THz respectively. A con- 
venient way to model such a complex structure is to 
solve the Maxwell's equations with appropriate boundary 
conditions using Finite Difference Time-Domain (FDTD) 
formulation ^,23^,27,]. A 3D-FDTD simulator with Uni- 
axial Perfectly Matched Layer (UPML)J23 surrounding 
the photomixer is developed for this purpose. Along the 
z-axis, the active volume consists of a few hundred nm of 
free space followed by DBRi , absorbing region with com- 
plex structure of normal metal electrodes, and pairs of 
low/high refractive index layers constituting DBR2 (Fig 
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FIG. 2: Optical electric field amplitude (a) and dc electric 
field strength (b) inside the D = 230 nm LT-GaAs absorbing 
layer of the photomixer design of Fig 1 optimized for Ao = 850 
nm. Incident optical field amplitude E^ = 1 V/m and T4 = 
±1 V. For clarity, fields in only the x — ±600 nm region along 
the y = plane are plotted. The grey rectangular regions 
in the plots show the cross-section of the interdigitated W 



lines with Lx 



7000, Ly 



2500, d = 10, 21 = 475, Z2 



555, z:i = 635, p = 300, w = 100, h = 1500, h = 1400 nm 
respectively. C = 4.90 £F and Rs — 2.4 Q for this device. 
The FP cavity is formed by three layers of (Ti02, SiaNx, 
CaF2) constituting DBRi (refractive indices (2.4, 2.0, 1.23 to 
1.36) [231 and thicknesses (100, 150, 190) nm for each layer 
respectively), four pairs of (Al20a;, GaAs) layers |24JI forming 
DBR2 (refractive indices (1.6, 3.53 — 0.068i) and thicknesses 
(130, 60) nm), and D « A absorbing LT-GaAs layer with 
refractive index 3.77 — 0.068J. 



la)|24^. For a given D, the refractive index and thick- 
ness of DBR layers are first optimized by calculating the 
reflection and transmission coefficients of the FP cavity 
using matrix methods|25|. The FP cavity is excited by 
a linearly polarized, plane wave propagating along +z 
direction, with gaussian (elliptic) intensity profile in the 
X — y plane. The source plane is placed in the free space 
above the DBRi. Frequency-domain Lorentz dispersion 
model|2y| is adapted for metals with negative dielectric 
constant such as W or Pt. 

Figure 2a shows the FDTD results of a three layer 
stack of interdigitated W lines embedded in 230 nm LT- 
GaAs absorbing layer of a photomixer design of Fig 1 
optimized for Aq = 850 nm. The plot clearly shows the 
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FIG. 3: Carrier transit time distribution ntr(ttr) for electrons 
(e) and holes (h) in the photomixer design of Fig 1. The 
parameters for red (solid) and blue (dot) curves are same as 
in Fig 2 with w = 100 nm where as for grey (dash-dot) curve 
only polarization is changed to Ey. For the green (dash) and 
cyan (short dash) curves {Lx, p, w) are changed to (7560, 
300, 60) and (7480, 320, 120) nm respectively while rest of 
the parameters remaining unaltered from those of Fig 2. The 
photomixer is excited by an 850 nm source at Pq — 60 mW for 
all structures. For clarity, the curves are shifted along y— axis 
and the relative shift can be easily obtained from the values 
at ttr ~ fs. 



FIG. 4: Calculated THz output power Pf from the pho- 
tomixer design of Fig 1 excited by Ao = 850 nm radiation 
with incident power Po = 60 mW. r = 4 ps for designs SI 
to S6 where as for S7 r = 1 ps. Parameters for SI are same 
as in Fig 2 whereas for S2 z^ — 645 nm, for S3 23 — 645 nm 
and incident radiation is Ey polarized, for S4 (p, w) are (320, 
120) nm respectively, for S5 {Lx, w) are (7560, 60) nm re- 
spectively, and for S6 {zi, 2:3) are (455, 655) nm respectively. 
The remaining unlisted parameters are the same in all cases. 
In all designs the interdigitated lines in different layers have 
no relative orientation angle differences. 



near field enhancement resulting from the thin normal 
metal electrodes. In the absence of the W electrodes, we 
obtain electric field amplitude maxima (maximum am- 
plitude 1.4 V/m) corresponding to three antinodcs in 
the standing wave formed between the DBRs. The 3D 
internal dc electric fields and C are computed by solving 
Laplace's equation with appropriate boundary conditions 
using FD techniques 123. Figure 2b shows the contour 
plot of the static electric field strength within the LT- 
GaAs layer. The data show that unlike the traditional 
MSM structure, the field strength in this design is well 
above the critical field (~5 kV/cm). Moreover, the field 
strength is ~90 kV/cm between neighboring electrodes, 
in particular at the center of the device where most of the 
carriers are generated due to gaussian intensity profile of 
the incident laser beams. In the rest of the volume it has 
a broad peak at 18 kV/cm (not shown). Results show 
that C — 4.90 fF and Rs = 2.4 fl based on the resist ivity 
of thin annealed or epitaxial W films (^5 /xf^— cm)J28j. 
The highest electric field inside the device is about four 
times lower than the breakdown field (500 kV/cm); there- 
fore device failure due to electric breakdown is unlikely 
for this Vfc = ±1 V. 

Photocurrent ip is calculated by first computing the 
electron/hole transit time distributions nl^(ttr), n%{ttr) 
shown in Fig 3 for the entire absorbing volume based 
on the above FDTD and FD calculations. From the 
data available in literature for GaAs and LT-GaAs, for 



ttr < tonset = 100 fs, the Carrier motion can be approxi- 
mated by ballistic transport with electron effective mass 
TOg — 0.088?TJo, where tuq is electron rest mass. This 
value is consistent with the slope of the linear portion 
of transient drift velocity curve obtained from Monte- 
Carlo calculations [l8J. The corresponding effective mass 
for holes is nih — O.IOGtoq. Effective masses larger than 
the accepted values for GaAs (0.063, 0.076 mo for elec- 
tron and light holes respectively) are considered so that 
longer ttr, thereby lower estimate of ip results from these 
calculations. In the photomixer of Fig 2, pure ballis- 
tic transport is applicable to electrons generated close 
to the -l-ve electrodes and holes generated close to -ve 
electrodes. These carriers generated predominantly in 
the near field enhancement region, transit through non- 
uniform dc fields in the 5 to 90 kV/cm range. It should be 
noted that tonset — 100 fs is considerably lower than the 
theoretical limit of ballistic motion in GaAs for this field 
strength range 16^. For ttr > tonset, electron motion is 
approximated by qasiballistic transport with time depen- 
dent velocity distribution similar to ref. 18 up to ttr = 3 
ps, and equilibrium drift velocity ('^ 1 x 10^ cm/s) for 
ttr > 3 ps. Hole motion is approximated by one third of 
the electron velocity at any given ttr > tonset resulting in 
velocities lower than those reported for GaAs|17J. 

The integral of curves in Fig 3 give the carriers {Ne, Nh 
for electrons and holes respectively) generated per period 
(T ~ 2.83 fs) of the 850 nm source with Pq = 60 mW. 
For the parameters of Fig 2, ^--^98% of incident power is 



absorbed in LT-GaAs layer while ritr has peaks at ttr ~ 
65 and 70 fs for electrons and holes respectively (solid 
red and dotted blue curves). Such a sharp peak followed 
by several satellite peaks and a long tail of the ritr distri- 
bution can be understood from the fact that most of the 
carriers are generated close to the electrodes in the near 
field enhancement regions where the static fields are also 
strong. The satellite peaks result from periodicity of the 
electrode structure and specific choice of the velocity dis- 
tribution. As ttr < 1 ps for majority of carriers in Fig 3, r 
is not a critical factor in determining the performance at 
THz frequencies. A sharp drop in ritr by a factor of 2 for 
electrons and a factor of 3 for holes at tj,. « tonset is due 
to the lower estimate of carrier velocities for ttr > tonset 
resultin g fro m the uniform field qasiballistic distribution 
function [13] applied to a case where fields are inherently 
inhomogeneous. Although, the above choice of velocity 
distribution for quasiballistic motion is consistent with 
the fact that over a large volume fraction of the absorb- 
ing layer the static field is ~18 kV/cm (not shown), a 
rigorous calculation should include quasiballistic velocity 
distribution appropriate for inhomogeneous fields. This 
would require ensemble Monte Carlo and carrier trajec- 
tory calculations carried out simultaneously. However, 
above calculations are adequate for obtaining the lower 
limit of the THz power output because a rigorous calcu- 
lation would probably produce a more uniform ntr{ttr) 
distribution around ttr ~ tonset without altering the dis- 
tribution for ttr < tonset Or the integral of ritr (ttr)- This 
will further reduce the number of carriers in the long tail 
of the rifr distribution. 

Based on above ntr{ttr) distributions, the number of 
electrons captured in the metal electrodes at time t will 
be 
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(t) 



<e„(*-)<r(i-ic)e-(*-*=)/"=dtc, (2) 
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where ri'i^^{t) — -^^^[1 4- cos(27r/i] is the number of elec- 
trons generated in LT-GaAs layer per second and tc is the 
carrier creation time. Similar expressions can be writ- 
ten for number of holes captured n'^api^) ^^id the elec- 
trons/holes available for conduction in the photomixer 
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e(»^cap(0 + 

Rl = 72 n 



i{t). In the calculation of ricap, riavi, ip{t) — 
(t)) (e is electron charge), and P/, we set 



cap 

T = 



2Te = 2r/j/3[lJ, and varied r in the 0.5 
to 6 ps range. For the parameters of Fig 2, the steady 
state electron density n obtained from the dc value of 
nl^i (t) for r = 4 ps is ~ 6 x lO^^ cm'^ (at Pq = 60 mW) 
with holes 1.5 times more numerous than electrons. We 
estimated that dipole fields arising from this space charge 
are negligible (< 2.5%) in comparison to strong fields (Fig 
2b) present in the photomixer. 

We have carried out calculations for over 
70 configurations by systematically varying 



the lasers {E^ or Ey), orientation of interdigitated lines 
in different layers, and DBR parameters. Fig 4 shows the 
Pf values obtained from Eq 1 and 2 for some of the con- 
figurations in the 0.1 to 10 THz range. The data shows 
that Pf oc /"^■''^ in the 0.5 to 6.5 THz range in contrast 
to /^^ roll-off of Pf for a conventional photomixer [lj|. 
A recent nip-nip photomixer concept is shown to have 
/-2 roll-off for / <1.5 THzfl^. Therefore, the design 
of Fig 1 exploiting SPPs offers significant improvement 
over the existing photomixer designs. Based on 3D 
FD computation of steady-state heat equation with 
appropriate thermal boundary conditions for the device 
parameters of Fig 2 (Pq = 60 mW, VJ, = ±1 V), internal 
temperature (T^) of the device is estimated to be ~200 
K above the substrate temperature requiring substantial 
cooling (to 77 K) when Pq is high. Thermal conductivity 
(fc) of various DBR layers and that of LT-GaAs are 
approximated by assigning k];'^''^''''^'' = 46 Wm^^K 

and fcV^-^"^^ 



lDBR 



I 

^DBR 



10 Wm-^K- 



where fc|| /k± are in-plane/out-of-plane conductivities. 
This internal heating therefore limits Pq to '--^60 mW. 
However, the structure of Fig lb, which is equivalent 
to two capacitors in series coupled to the antenna or 
transmission line, offers an alternative. If it is excited 
by TEMoi mode lasers with total power Pq = 120 mW 
as shown, the output can be increased to 2P/ (of Fig 
4) without worsening Ti or n. To demonstrate the 
applicability of this method to other semiconductor 
materials, similar calculations are carried out for two 
more photomixers based on Be doped Ini_a:Ga3;As and 
GaN optimized for operation at Aq = 1550 and 343 
nm respectively. Both structures show strong near 
field enhancement arising from SPPs similar to Fig 2a 
and further work is underway to calculate the device 
efficiencies. Author wishes to thank John Goodkind for 
introducing him to the fascinating subject of Auston 
switches. 
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